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Abstract 

In this work we study the nonlinear dynamics of the static and the driven ellipse. In the static 
case, we find numerically an asymptotical algebraic decay for the escape of an ensemble of non- 
interacting particles through a small hole due to the integrable structure of the phase space of the 
system. Furthermore, for a certain hole position a saturation value in the decay that can be tuned 
arbitrarily by varying the eccentricity of the ellipse is observed and explained. When applying 
harmonic boundary oscillations this saturation value caused by librator type orbits is gradually 
destroyed via two fundamental processes which are discussed in detail. As a result, an amplitude 
dependent emission rate is obtained in the long time behavior of the decay, suggesting that the 
driven elliptical billiard can be used as a controllable source of particles. 
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I. INTRODUCTION 



Billiards belong to the most widely studied Hamiltonian systems. They possess many 
classical and quantum mechanical properties of complex dynamical systems if, [2], [J. More- 
over, models of statistical mechanics can be reduced to billiards 4]. For example, one of 

the simplest billiards, particles inside a rectangular box, is an idealization of the physical 

Tl 

situation of nucleons confined inside a nucleus [5|. Mathematically rigorous studies of bil- 
liards go back to the early seventies, e.g. Bunimovich proved that stadia are ergodic Q], 
using concepts developed by Sinai. In recent years, a renewed interest in billiards has come 
up, due to the possibilities of realizing them experimentally, for example by using ultracold 



atoms confined in a laser potential [71], microwave billiards 



or mesoscopic quantum 



dots [10]. Even for the design of directional micro-lasers, billiards are relevant 



11]. Besides 



this, interesting theoretical results were obtained, including a justification of a probabilistic 



approach to statistical mechanics 



12j. Very recently, it has been shown [13j, that a con- 



nection exists between billiards and one of the major unsolved problems in mathematics, 
the Riemann hypothesis: the authors found an analytic expression for the escape rate of a 
circular billiard with two holes, involving a sum over the zeros of the Riemann zeta function. 

A natural generalization of billiards with a static boundary is to apply a driving law 
to the billiard wall. For instance, 'Bohr's liquid drop model' from nuclear physics can 
be regarded as a time-dependent billiard [ijj]. For this simple looking model still many 
questions remain open 151 ] . Another example is in plasma physics, where time-dependent 
billiards represent models for acceleration of particles in a magnetic bottle, see and ref. 
therein. In conclusion, there are many branches of physics in which billiards, specifically 
time-dependent billiards, serve as models for more complex systems, capturing the key 
features and behavior of the original problem. 

Ultracold atoms in a billiard formed by beams of light allow for the possibility of gener- 
ating arbitrary geometries and changing them in time, as well as varying parameters such 
as beam width, softness of the potential etc. in time. Of special interest is the possibility to 
probe the dynamics by analyzing the escape rates [l?!, Q, Q, 3; which has up to date 
only been performed for static billiards. Introducing noise and decoherence and studying 
the role of quantum and many-body effects are further intriguing goals [7]. 

Regarding time-dependent billiards, there exist several investigations in the literature 
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281 ] . A crucial question for these systems is whether Fermi 



acceleration occurs or not. This is examined in refs. [211 . [22], [23|, [24| and very recently in ref. 
[27I . In 2l| it is shown that when using smooth forcing functions, the existence of invariant 
spanning KAM curves in phase space limits the energy gain of the particles, whereas non- 
smooth forcing functions, especially random oscillation lead to unbounded energy gain, see 
also ref. [27J and references therein. In ref. 22], the authors conclude with the hypothesis: 



"A random element in a billiard with a fixed boundary is a sufficient condition for the Fermi 
acceleration in the system when a boundary perturbation is introduced. " . 

Within the existing studies of classical time-dependent billiards only little emphasis is put 
on systems with a finite horizon and (to our knowledge) none to the corresponding escape 
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26|]. In ref. [l6j, the 



rates. Very few works deal with the time-dependent ellipse 
average velocity as a function of time and the Poincare surface of section of the dynamics of 
the ellipse for different driving laws are studied numerically. Depending on the driving laws 
and the initial velocity of an ensemble, the integrable structure of the phase space is more 
or less destroyed compared to the static case. The velocity of the particles stays bounded 

'he authors point out, that all conditions 



291 ] which predicts this boundedness of the 



in all cases, i.e. no Fermi acceleration occurs, 
are satisfied in order to apply Douady's theorem 
velocity. A mathematical study of periodically driven ellipses is given in |25j. The authors 
show that in principle it is possible to destroy the diametral 2-periodic orbit via boundary 
oscillations and give strong evidence that the opposite - stabilizing an unstable periodic 
orbit with the use of driving - is not possible. 

The above discussion shows that essentially little is known about escape rates in classical 
time-dependent billiards. Apart from being of fundamental interest this type of driven 
dynamical systems is nowadays well within the reach of experiments, as indicated above. 
Moreover our investigation will demonstrate that time- dependent billiards might provide 
us with a tunable source of particles. As we shall see, the escape rate and the velocity 
distribution of the escaping particles strongly depends on the driving properties, such as 
amplitude and frequency of a periodic driving. 



In this work, which is an extension of our recently published letter [30j, we focus on the 
driven elliptical billiard. Its static counterpart is integrable, due to the existence of a second 
constant of motion, the product of the angular momenta around the foci. Thus the phase 
space of the ellipse possesses a more complex structure, consisting of librators and rotators, 
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than the prototype integrable billiard, the circle. Naturally, it would be also interesting 
to study driven billiards whose static counterparts have mixed or chaotic dynamics. Yet 
the clear partition of the phase space into librators and rotators of the ellipse simplifies the 
analysis in the presence of the driving considerably. This allows us to study e.g. transition 
between librator and rotator orbits and to discuss associated physical phenomena in an 
intuitive way. 

This article is structured as follows: In section [Til we discuss fundamental properties and 
escape rates in the static ellipse. The generalization to time-dependent ellipses is treated in 
section IIHI The fundamental processes leading to the destruction of the librator orbits are 
displayed in section [III CI followed by an analysis of the angular momentum, section ITV] and 
the velocity, section [V] Finally, a summary is given in the last section. 

II. STATIC ELLIPSE 

A. Fundamental Properties of the Dynamics in the Ellipse 

In a two-dimensional static billiard, the orbit of a particle can be completely specified 
by providing the sequence of its positions Sj (measured by the arclength) or (ft (see eq. [I]) 
on the boundary B and the directions Pi = cosctj immediately after each collision, since 
the particles travel ballistically in between collisions, where <X{ is the angle between the 
forward pointing tangent and the velocity of the particle at the i-th collision point. The 
corresponding discrete mapping A4 is area preserving in the phase space variables s and p 
The boundary B of an ellipse is given by 



with A > B > 0, thus A and B being the long and the short half-diameter respectively. The 



In anticipation of the time-dependent problem, we describe the direction of a particle by 
its velocity v = (v x ,v y ). If we demand without loss of generality \v\ = 1, there is a one to 
one correspondence between the velocity v and p at the collision points. At a certain time t, 
the position of the particle starting at t = at x = (x , y ) G B with the velocity (v X fi, v Vt o) 



B = {x(f) = A cos ip, y((p) = Bsm<p)\0 < ip < 2n} 



(1) 



dimensionless numerical eccentricity can be written as e 



y/l-B 2 /A 2 . 
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is given by 



x{t) 

y(t) 



Xo + v xfi t 

VO + Vyfit. 



The particle will hit the boundary at X\ at the time t±. 

2B 2 x v xfi + 2A 2 y Q v yfi 



(Av yfi )> + (Bv xfi ) 2 



(2a) 
(2b) 

(3) 



To get the new velocity V\, we parameterize xx e B by ipx and calculate the inward pointing 
normal vector ni, |ni| = 1 at (fx. This results in 

vi = v -2(n a -v )-ni. (4) 

Equation (jl]) can be easily extended to time-dependent boundaries, see section IHIt where 
momentum transfer from the moving wall to the particle takes place. 

1 




FIG. 1: PSS of the ellipse, (upper part) and typical trajectories (lower part), A = 2, B = 1. The 
rotator orbit repeatedly touches a confocal ellipse, the librator orbit a confocal hyperbola. 

The dynamics in the ellipse is completely integrable, see Fig. [TJ In addition to the 
energy, there is another constant of motion F(<p,p), restricting the orbits to invariant curves 
in phase space: 
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l + (l-£ 2 )cotV~£ 2 

F(<p,p) can be interpreted as the product of the angular momenta (PAM) about the two 
focus points [1]. There are two different types of orbits, rotators and librators, in the ellipse 
separated by the separatrix, see Fig. [IJ Librators cross the x-axis between the two focus 
points and touch repeatedly a confocal hyperbola. In the PSS, they appear as deformed 
circles around elliptic fixed points, exploring a limited range in p as well as in (p. Rotator 
orbits travel around the ellipse, exploring every value of p, but only a small range in p 
(except if they are very close to the separatrix), repeatedly touching a confocal ellipse. 

In terms of F(p, <p), we can distinguish tree different cases: 

1. F(ip,p) > corresponds to the rotator orbits with an elliptical caustic. 

2. F(<p,p) < are the librators with hyperbolic caustic. This includes the two elliptic 
fixed points at (cp = tc/2, p — 0) and (ip = 3ti/2, p = 0) corresponding to a period two 
orbit along the minor axis with F(7r/2, 0) = — e 2 / (1 — e 2 ). 

3. F(<p,p) = corresponds to the period two orbit along the major axis with (p = 0, p = 
0) and (<p = 7r, p = 0), seen in the PSS as two hyperbolic fixed points. 

The topology of the PSS, Fig. [IJ is dominated by two isolated periodic orbits. The 
condition for stability is according to ref. l| 

n > unstable 

1 (6) 



2R(p) 



< stable, 



where R(p) is the radius of curvature and p is the distance in coordinate space between two 
successive collisions. Using this stability criterion, we get for the periodic orbit along the 
long diameter p/2R = 1/(1 — e 2 ) > 1 and therefore it is unstable. In contrast, the orbit 
along the short diameter obeys p/2R = 1 — e 2 < 1 and is stable. 

n 

According to ref. [16], Poncelet's theorem on projective geometry can be applied to 



elliptical billiards [3jJ. It states, that all trajectories possessing the same value of F(p,p), 
share the same caustic and the same type of dynamics. In the case of periodic orbits this 
means, that given one periodic orbit with a certain value of F((p,p), every trajectory with 
the same value of F(<p,p) is also periodic and has actually the same period. Consequently, 



the only isolated periodic orbits are the two discussed two-periodic orbits, all the other 
periodic orbits are non-isolated and form families. 

B. Escape rates 

Let us focus now on the escape rates of a static elliptical billiard with a hole placed 
on its boundary. In this subsection we use for all simulations A = 2, B = 1, i.e. the 
numerical eccentricity e = \/3/2 « 0.87. The number Nq of particles in the initial ensemble 
is iVo = 10 7 . Each particle is propagated at most 10 6 boundary collisions unless it does not 
escape earlier. The initial conditions (ipo, «o)j, « = 1,2,... 10 7 (the index i stands for the ith 
particle) are chosen randomly. Note that the angle «o is distributed uniformly in [0, n], not 
p = cosa . We choose two different hole positions <p& = and <p& = n/2. The hole size A 
is set to A = 0.03 (measured in (p). (p& = corresponds to a hole lying in the very right of 
the ellipse of Fig. Q], and v?a = tt/2 corresponds to the location at the very top of the ellipse. 
The reason for this choice is the following: If the hole lies at (p& = 0, none of the librator 
orbits can escape, since their invariant curves are not connected with the hole, whereas if 
V^a = 7r /2, all orbits can participate in the decay. In both cases, all rotator orbits can escape 
(as long as they are not periodic), since they are ergodic with respect to the phase space 
variable ip. The main data of these simulations is the number of remaining particles in the 
billiard as a function of the number of collisions N(n) or the elapsed time N(t). Note that 
we refer to N(t) as the escape rate, as done in the literature, whereas we will call N(t) the 
emission rate. 

The results of the simulations are shown in Fig. [2j Two qualitatively different behaviors 
of the decay are observed for the two different hole positions: 

1. If sa = 7r/2, the fraction of remaining particles N(n) as a function of the number of 
collisions approaches zero for n — > oo. 

2. If sa = 0, N(n) approaches a saturation value N s (e) > after roughly n — 2 ■ 10 3 
collisions. 

The saturation value N s (e) in the case sa = is of course caused by particles traveling on 
librators. Since the librator orbits are not connected with the hole, these particles will stay 
forever in the billiard. We will derive an exact expression for N s (e) in the next section. In 
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FIG. 2: Semi-logarithmic plot of the escape rate, two different hole positions are shown, double- 
logarithmic scale in the inset. 

both cases of the hole position, the short time behavior, Fig. El of the decay is exponential 
N(n) ~ exp(— rn) (roughly for the first 50 collisions in the case sa = and 300 collisions in 
the case sa = 7r/2). The decay constant r is approximately given by r w A/27T 32]. The 
long time behavior (n > 2 • 10 3 ) of N(n) in the case sa = 7t/2 corresponds to an algebraic 
decay N(n) ~ n~ c , seen as a straight line in the inset of Fig. El This power law decay is 



typical for integrable systems and known in the literature, see e.g. ref. [5j] or 



but there 



is no work discussing the case of the ellipse, except for ref. 19j, where the algebraic decay, 
even though not in such detail, is observed experimentally. 

A heuristic model explaining this algebraic behavior is provided in The discussion 
given there holds for a rectangular box, where |p • e n | (e n is the unit vector normal to the 
opening) is a constant of motion. Nevertheless, the results obtained there can be easily 
transferred to the case of the ellipse by replacing |p • e n | by F(ip,p). According to [5], the 
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fraction of remaining particles should decay for large n like N(n) ~ n~ l . The extracted 
value from our data is N(n) ~ n~ im for n > 3 ■ 10 3 , i.e. in very good agreement with the 
analytical prediction. 



C. Saturation value N s (e) 

Let us now study whether the escape rates depend on the numerical eccentricity e. Indeed, 
the qualitative behavior of the decay remains unchanged, only the saturation value N s (e) is 
different for different values of e for the hole at the short side of the ellipse, sa = 0. This 
becomes immediately clear if one considers that e determines the degree of deformation 
compared to the circle: Since in the circle there are exclusively rotator orbits, N s (0) should 
be zero and with increasing e the offset N s (e) should increase too. 




Numerical eccentricity £ 

FIG. 3: Dependence of the saturation value N s on e. 
The fraction of escaped particles as a function of e is shown in Fig. [31 One can also cal- 
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culate N s (e) from the following theoretical considerations, the result is excellent agreement 
with the numerical data points, see Fig. [31 

All initial conditions corresponding to particles propagating on librators lie inside the 
area A n (e) bounded by the separatrix of Fig. [H We denote by f(e) = A n (e) / A PS s the 
ratio of Ajj(e) and the total area Apss of the phase space. 

ApSS = (Pmax — Prnin) ' (<Pmax ~ <Pmin) = 4vT. (7) 

To calculate A/j(e), we need an analytic expression of the curve belonging to the upper half 
of the separatrix, that is a function p sx (<p). Then, this area is given by 

A n (e) = 2 / d<pp sx (<p). (8) 
Jo 

We know that for the motion along the separatrix F(<p,p) = 0. Thus, we can exploit (jSJ) 
and get 

P-M=^ = 0) = J f (9) 



We see immediately that Ajj{e) (jSJ) depends on e, and so does the saturation value. To 
obtain N s (e), we have to account for the fact that the initial conditions are distributed 
uniformly in the a, tp - and not the p, tp - space. Hence, 



e 2 



a sx ((p) = arccos(p sx ((p)) = arccos W - — -5— (10) 

" 1 + (1 — e z ) cot (p 

and as a result f'(e) = A' II (e)/2ii 2 , similar to ([7]) and A' n (e) is 

A' n {e)=2 d<pa sx (<p). (11) 
Jo 

The fraction of escaped particles is just 1 — /' and the saturation number is 

N s (e) = f(s). (12) 

In Fig. [31 perfect agreement between the above presented analytical considerations and the 
numerical simulations can be seen. As a consequence, varying e allows us to control the 
number of particles being emitted. 
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III. TIME-DEPENDENT ELLIPSE 



In this section, we investigate the escape rates for the time-dependent ellipse. Since the 
boundary transfers momentum to the particles upon collisions, their energy is not conserved 
any more. The collision point of a particle with the boundary is not defined by ip only, but 
we need additionally the time t to make the point well-defined in coordinate space, since 
the boundary B(t) depends explicitly on t. Likewise, the direction of a particle has to be 
described by v = (v x ,v y ) and not just by p = cos a, since \v\ ^ const. Representative 
visualizations of the resulting 4D-phase space, like the 2D-PSS for the static billiard, are 
difficult to achieve, if not to say prohibitive. 

To drive the ellipse, we apply harmonic oscillations to its boundary B(t) 

B(t)={b(ifi,t)\^e[0,27r)} (13) 

b^t)=( X ^ t) ) = ( A{t)COS A (14) 
\y(<p,t) J \B(t)sm<p J 

where A(t) and B(t) are given by 

A(t) = A + C sm(ut + 5) (15a) 
B(t) = B + C sm(ut + 5) (15b) 

C > is the driving amplitude and 5 is a phase shift. Ao, Bq and C have to be chosen in 
a way, that A(t) > and Bit) > for all t. We refer to ( floT) as the breathing ellipse. As 
already done in section [Til we set A = 2 and B = 1, and use values of C between 0.01 and 
0.30 only. The velocity u(ip,t) of the boundary and the numerical eccentricity are 

uCcos(ujt+5)cos(p 
u[<f,t) =| I , (16) 

ujC cos(u;t + 6) sin ip 



£(t) -^ 1 "(2 + Csin(o;t + W 



A. Mapping 



Just like in the static case, a discrete mapping is sufficient to characterize the full dynamics 
of a particle. Consequently, the trajectory of a particle consisting of N bounces is given by 

C = {(t , <p Q , «o), (h, (fl, Wi), . . . , (t N , V?7V, v N )}. (18) 
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The mapping for the next collision time t n+ \ is determined implicitly by 



v™{t n+1 -t n ) + x n \ 2 ( v"(t n+ i - t n ) + y, 



2 



1 = (19) 



A + C sm(ujt n+1 + 8) J \B + C sin(o;t n+ i + 8) 

where for a given t n and (p n , x n and y n are calculated from (jHj) and t n+ i is defined by 
the smallest t n+ i > t n that solves (fT9l . The next collision point is given by x n+ \ = x n + 
v n (tn+i — t n ) and <p n +i can be obtained by inverting (|T4]) . Once (t n+ i,<p n+ i) is determined, 
the next velocity v n+ i is given by 

v n +i = v n -2 [n„ + i • (w„ - u n+1 )] • n n+ i (20) 

where the boundary velocity n„ +1 is given by (flEj) and the the normal vector by n n+i = 
<+i/K+il> <+i = (-B(t n+1 )cos(f n+1 ,~A(t n+1 )sm(f n+1 ) T . 

The maximal velocity change of a particle upon a single collision with the boundary is 
according to eq. ( |20l) A|«| = ±2co>C. Thus, it can happen that the particle undergoing 
a boundary collision at the time t' is not reflected, in a sense that the sign of the velocity 
component normal to the boundaries tangent is not reversed, but continues traveling outside 
B(t'), but of course still inside B(t > t'). This is the case if the ellipse is expanding and 
u n < v n < 2u n holds, where v n and u n are the normal component of the particle and the 
boundary velocity before the collision. As a consequence, the angle a between the tangent 
t and the velocity v is not restricted to the interval [0, n] as it was in the case of the 
static ellipse, but now a G [— n, ir]. Upon such collisions with the expanding boundary, the 



particles are always slowed down, they lose energy 3J, |35(, whereas upon collisions with the 
contracting ellipse they gain energy. 



B. Escape Rates 

We focus on the case sa = 0, in order to examine the effect of the driving on the number 
of particles in the billiard. For the static case, the saturation value was caused by the librator 
orbits, see section III CI We assume that these librator orbits will be deformed or partially 
destroyed by the driving, leading to a non-vanishing decay even for large times. On the 
other hand, we expect no stabilization of the rotator orbits, i.e. no deformation in a way 
that they will not escape. All periodic orbits become unstable when applying the driving 
and in ref. [25| it was concluded, that it is impossible to trap unstable periodic orbits in the 
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ellipse via boundary oscillations. Note that this is not true in general for driven systems, 
unstable periodic orbits can be stabilized by a driving force, e.g. in the Kapitsa pendulum 



36]. 
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FIG. 4: Fraction of remaining particles in the IVE as a function of time for different values of the 
amplitude C, semi-logarithmic plot in the inset. 

We consider two different borderline cases: \v \ « ujC (intermediate velocity ensemble 
(IVE)) and \vq\ ^> ujC (high velocity ensemble (HVE)). In the first case, the velocity has 
the same order of magnitude as the the boundary velocity. This leads to a momentum 
transfer (maximal 2ujC) of the same order of magnitude compared to the initial momentum 
and we expect significant changes in the dynamics. In the second case, the particles move 
much faster than the boundary, consequently the momentum transfer will be very small, the 
dynamics will be similar to the one of section [III Naturally, it would be also interesting to 
examine the case |«o| ujC. However, the first few collisions then accelerate the particles to 
velocities \v\ ~ ujC and, after a short time, we encounter the situation of the first case. The 
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FIG. 5: Same as Fig. [1] for the HVE, additionally 27r-oscillations of the decay are shown in the 
inset. 

parameters of the simulations are iVo = 10 5 , u = 1, C = 0.01, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30 
and 5 = 0. To ensure that all the particles move inside the billiard, we let them start 
on the smallest ellipse (for a given C), the initial position (po and the initial angle ao 
are chosen randomly. The initial velocity Vq is given by vq — (cosao,sino;o) (IVE) and 
vq = 100 • (cosao, sinao) (HVE) respectively. The fraction of remaining particles Nc(t) as 
a function of time for different amplitudes C is shown in Figs. H] (IVE) andO (HVE). 

Firstly, we describe the behavior of the IVE. We observe a short but fast decay (t < 500), 
followed by a transient (500 < t < 5000) in which the decay slows down and for t > 5000, 
the decay is much slower than the initial fast decay for all values of C. At t = 10 4 , the 
values of the fraction of remaining particles are ordered according to the driving amplitudes, 
the lower C is, the higher is Nc(t = 10 4 ), i.e. Nc(t = 10 4 ) depends monotonically on the 
amplitude C. For t > 10 4 , Nc(t) does not stay constant, but is still decreasing. The absolute 
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value of the emission rate Nc(t) is the larger the larger C is. This can be seen nicely in 
the double- logarithmic plot of Fig. HI For values of t between 10 4 and 10 5 , we encounter 
approximately an algebraic decay Nc(t) ~ t~ w (we remark that this algebraic decay has 
been numerically shown to exist for much longer times than illustrated in Fig. Hj), where the 
decay constant w increases monotonically with increasing C (this fact is based not only on 
the four values of the driving amplitude C shown here, but on simulations carried out for 
20 values of C between 0.01 and 0.30). 

The subdivision of the behavior, into fast initial decay - transition period - slow (near 
algebraic) decay, is even more pronounced in the case of the HVE, see the inset of Fig. |5j 
An exponential decay for small values of t slows down at around t ~ 5 and the fraction 
of remaining particles seems to approach a constant value. From Fig. [5] however we see, 
that the fraction of remaining particles decays for t > 10 roughly according to an algebraic 
decay (at least for small values of the driving amplitude) N(t) ~ t w with a decay constant 
w. If we compare the fraction of remaining particles at t = 50 for different values of the 
amplitude (inset of Fig. [5]), we see that they are monotonically ordered according to the 
driving amplitudes. Surprisingly, most of the particles remain within the billiard in case 
of the largest driving amplitude C = 0.30 and the smallest fraction remains in case of the 
smallest amplitude C = 0.05. The explanation of this effect is provided later, in section HVl 
when we examine the dependence of the PAM F(ip,p) on the driving. 

In the inset of Fig. [51 a modulation of the escape rate with period T = 2tt can be 
seen, being exactly the period of the applied driving law ( 1T4"1) . Specifically, for t > 10, 
where all particles starting on rotator orbits have already escaped, Nc(t) ~ const, during 
approximately 11/12 (empirically observed) of one period and subsequently Nc{t) 7^ 
during a time interval T/12 only. From this behavior it is evident that the ellipse operates 
from a certain time on as a pulsed source of particles. These repeated intervals are centered 
around points t m of maximal extension of the ellipse, t m = (4m + l)vr/2, m = 2,3,4, .. . 
During the expansion period, dominantly vertical but also horizontal processes turn librators 
into rotators. The moving ellipse remains for a comparatively long time period in the vicinity 
of the extremal configuration at t m and consequently the newly created rotators escape. 
Therefore, the dynamics is effectively probed during these short time intervals centered 
around t m . During the contraction period, the librators are stabilized via vertical processes, 
consequently Nc(t) ~ during 11/12 of a period T. 
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C. Mechanisms for the destruction of the Librators 

In the driven ellipse librators can escape from the billiard whereas this is not the case for 
the static ellipse. There are two fundamental processes that perturb or even completely de- 
stroy the librator orbits (unprimed variables denote the static, whereas primed ones describe 
the driven system): 

1. Vertical process: The angle of incidence of a collision does not coincide with the 
reflection angle because of a change of momentum due to the motion of the boundary 
of the ellipse. In phase space, the momentum then undergoes a certain change Ap 
upon a collision and the particle moves vertically in the PSS. 

2. Horizontal process: A particle that would hit the boundary at a certain point ip, hits 
the boundary in the driven case at ip', simply because the ellipses boundary has moved, 
whereas p stays nearly unchanged. This corresponds to a horizontal move in the PSS. 

These processes are fundamental in the sense that every change AF can be decomposed 
(at least for small changes (A(p, Ap)) into AF = AFh + AF V , where AF^„ denote the 
individual changes caused by the horizontal, the vertical process, respectively. 

In general, these effects do not appear isolated but a combination (A(p, Ap) of both 
will occur in a single collision. We can compare the orbits p9 of the driven ellipse to 
the corresponding ones of the static ellipse (<Pi,Pi) by considering the quantity F(ip,p) (see 
eq. <^). In contrast to the case of the static ellipse where F(ipi,Pi) = const. Vz we have 
F{<p'i,p'j) 7^ F{y>'jiP'j) (i 7^ j) f° r the driven case, i.e. F is no longer a constant of motion. 
The difference AF (see Fig. [1]) upon a collision is a measure of whether a librator approaches 
the separatrix (AF > 0) or whether it moves in phase space towards the position of the 
elliptic fixed points (AF < 0) of the static case. An increase with respect to F reflects the 
dependency of 'moving' in phase space from confined librator to escaping rotator orbits. 

1. Vertical Processes 

To isolate this effect, we examine a particle that hits the boundary at ip = tt/2 under 
a certain angle a in the static case. The velocity of the particle can be written as Vq = 
{—v cos a, v sin a), v = \vq\. In the driven ellipse, we will assume that the particle hits the 
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FIG. 6: Vertical process in coordinate space. 

boundary in the neutral position (A(t) = A , B{t) = B ), so we have ip' — if — ir/2, =>- 
Aip = 0. The boundary velocity u(ip,t) of the ellipse is maximal at this configuration 
and has a vertical component only, u n = u y = ±uC, depending on whether the ellipse is 
expanding "+" or contracting "-" . The of the particle at the next collision in the static case 
is Vi = (— cos a, — sin a) T and hence pi = Po = cos a. The corresponding velocity in the 
driven case is v[ = (—vcosa, — v sin a ± 2uC) T . Now, p' = p Q ^ p[, since 

. cos a: 



Pi 



AluC Aco 2 C 2 
1 =F sin a H — 



(21) 



/ 



If the ellipse is contracting ("+" sign in the factor /), p[ = cos a/ f is smaller thanpi = cos a 
because / > 1, i.e. Ap < (p[ — pi + Ap). In phase space the particles 'moves' towards the 
elliptic fixed points, which does not lead to the destruction of the librators. If the ellipse is 
expanding ("-" sign in the factor /), Ap will be larger than zero if / < 1. This is equivalent 
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and since v n = v sin a and u n = ojC eq. (J22l) is equivalent to i> n > u n , which is a necessary- 
condition for a collision to take place. If additionally to eq. (|22p sin a < 2uC/v, then 
v n = v y > and thus a' < 0, see Fig. 0, the particle is not reflected in the sense that the 
sign of v n is not reversed. For sin a > 2uC/v it follows that v n = v y < 0, the particle is 
reflected. In both cases Ap > 0, the particle moves towards or even beyond the separatrix; 
in the latter case, the librator has changed into a rotator. This process can therefore lead to 
the destruction of the librator orbits. In Figs. [6] and El the process is shown in coordinate, 
as well as in phase space. This process happens for any value of ip, nevertheless, the case 
cp ~ 7r/2 is for two reasons especially important: 

1. The absolute value of the change Ap is largest for cos a m sin a, for particles on 
librators this is approximately true if (p « 7r/2. 
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2. For a constant Ap, the corresponding AF(ip,p) is largest for ip = ir/2, because the 
vertical spacing of the invariant curves is smallest at (p — 7r/2. 

Consequently, the horizontal processes that contribute the most to changing a librator into 
a rotator occur mainly around (p ~ 7r/2 (and « 37r/2 because of symmetry). 

2. Horizontal Processes 
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FIG. 8: Horizontal process in coordinate space. 



We consider a particle that starts moving in a static ellipse corresponding to the smallest 
possible ellipse (given a certain amplitude) in the driven case, hitting the boundary of this 
static ellipse at ip = n/2 under a certain angle a, see Fig. El We choose all the parameters 
such that the particle will hit the boundary of the driven ellipse at the time t' at the 
position ip' ^ (p, when the boundary B(t') has its maximal extension. The new angle a' is 
approximately equal to a, see Fig. |HJ This rough estimate becomes better with increasing 
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FIG. 9: Horizontal process in phase space. 

a and decreasing distance between the two ellipses, thus Ap pa 0, i.e. the particle moves 
horizontally in the PSS, see Fig. [9j To calculate (p' or Ap as a function of C and a is 
very tedious and the exact result is not very helpful. We therefore restrict ourselves to an 
approximation and linearize the boundary of the ellipse locally at tp = n/2. The collision 
point on B(t) is 

. y (2 + C) ■ tana 
ip = arctan — pa arctan — . (23) 

In general we have to respect the sign of x and y to obtain ip' from (|23|) . |A<^| = \tp' — <p\ 
depends sensitively on a and decreases with increasing a. \A<p\ is largest for <p pa tt/2, since 
a there reaches its minimal value for the librator orbits. The sign of the corresponding 
change AF depends on the sign of Aip and the quadrant in which ip lies (e.g. for A(p < 
AF > if ip lies in the second or forth quadrant). There is no obvious region with respect 
to (p where the destruction of the librators occurs mainly, since two opposite effects balance 
each other: 
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1. The absolute value of the change \Aip\ is largest for ip »s 7r/2. 

2. Given a certain Ayj, the corresponding changes AF increase with increasing distance 
between (p and tt/2, because the horizontal spacings of the invariant curves is smallest 
there, see Fig. 

Nevertheless, the position (p = tt/2 (and <p = 3n/2) is in a way exceptional: any change 
Aip, independent of the sign, results in a change AF > 0, leading to a destruction of the 
librators. 



D. Qualitative Model of the Decay 

We isolated two processes that are able to change F(<p,p) in the course of the dynamics. 
Particles on librator orbits are scattered upon boundary collisions either towards the elliptic 
fixed points or towards or even beyond the separatrix. For a single particle, such a scattering 
process happens at every collision, the effective change AF after a certain time depends on 
the sequence of these processes, hence AF = AFi + AF2 + • • • + AF n after n collisions. 
This effective change in F(<p,p) is very difficult, if not impossible, to predict, since each 
individual change AFj depends on four parameters already: 1. the absolute value \vi\ of 
the particle velocity, 2. the angle otj of the velocity with the boundary, 3. the location (pi 
of the collision point on the boundary and 4. the time which determines the position of 
the boundary and the boundary velocity ttj(ii) of the ellipse (these four parameters are of 
course just the variables of the four-dimensional discrete mapping, see section [III Aj) . Now 



we consider not only a single particle, but an ensemble of N particles with initial conditions 
(ipj, Vj), j = 1,2,...N. The effective change (AF)j (where the index j indicates the jth 
particle) after n collisions can vary significantly from particle to particle, since the sequence 
of these four parameters will be very differently for each individual particles. Each of the TV 
sequences is governed by applying the discrete mapping of section HH Al n-times on each initial 
condition (ipj, Vj). The underlying nonlinear dynamics of this discrete mapping and the fact 
that all particles start from different initial conditions leads to such unique sequences and will 
cause consequently large fluctuations in the effective AF and accordingly large fluctuations 
in quantities that depend on F(<p,p). In the following, a qualitative explanation of the 
escape rates Nc(t) (decay) of the IVE and the HVE is given. 
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We focus on the HVE first. The initial fast decay of the number of particles (t < 5) is due 
to the rotator orbits that are connected with the hole and escape very rapidly. Additionally 
some of the particles starting on librator orbits near the separatrix F J$ contribute. 
The longer-time decay (t > 10) is caused by particles starting on librators that have been 
scattered across the separatrix. The closer an orbit of a particle lies near the elliptic fixed 
points, the longer it takes until the effective change AF is big enough to reach the separatrix 
(F = 0). From equations (1211) and (123]) it follows that the individual changes AFj under a 
single collision increase with increasing amplitude C . This explains the increasing emission 
rate Nc(t) with increasing C, since at a given time t, the number of particles that can 
participate in the decay is larger for larger values of C. The decay in the transient region 
(5 < t < 10) is caused by a superposition of the tail of the initial fast decay (roughly 
exponential) and the onset of the slow (roughly algebraic) decay. 

With very similar arguments, the decay of the IVE can be explained qualitatively. Since 
the velocity of the particles and the velocity of the boundary are of the same order of 
magnitude, the changes AF are much larger compared to the ones of the HVE. This leads 
to a very early onset of the slow (algebraic) decay, consequently the transient region is 
broadened. 



IV. ANGULAR MOMENTUM 

To validate the qualitative model of section [Til Dl we investigate the PAM F(ip,p) further. 
The contours of F((p,p) are shown in Fig. [TJ depending on the initial value of F, particles 
move on rotator or librator orbits, see section III A[ 

Through out the following sections, we analyze properties such as the escape time for an 
ensemble of particles possessing certain initial distributions in e.g. phase space or the PAM 
F. We emphasize that in case of the distribution of the F-values we always refer to the 



initial distributions of F at t — 



A. Escape Time versus Initial Conditions in Phase Space 

Exemplarily, the escape time as a function of the starting points in phase space is shown 
in Fig. HU1 (C = 0.10, IVE), i.e. we assign to each initial condition (<p , p ) (10 5 particles) an 



22 



x 10 




FIG. 10: Escape time versus initial conditions in phase space for C = 0.10 (IVE). The escape time 
is largest for initial conditions close to the elliptic fixed points. 



escape time. Large values of the escape time correspond to initial conditions belonging to 
librator orbits lying around the elliptic fixed points at (if = tt/2, p = 0) and (if = 3tt/2, p = 
0). On the other hand, initial conditions corresponding to small values of the escape time 
lie around areas which correspond to the rotator orbits. Overall, the results of Fig. [TU] are 
in good agreement with our predictions of section [III D\ where we derived large escape times 
for particles with initial conditions close to the two elliptic fixed points and short escape 
times for particles starting on librator orbits. 



B. Escape Time versus Initial Angular Momentum 

In this section, we investigate the escape time t esc versus the initial angular momentum 
F(if ,p ) of the corresponding ensemble of particles in phase space for different amplitudes 
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C. We consider escaped particles only. The results of the IVE are shown in Fig. [Hj 
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FIG. 11: Escape time versus initial angular momentum F(tpo,po) (F = corresponds to the 
separatrix in the static ellipse and F = —3.1 to the elliptic fixed points) in the IVE for C = 0.10 
and C = 0.01 (inset). 

For small escape times t esc < 10, there is a narrow, serpentine chain in which all the pairs 
(t esc , F) lie. Only pairs with F > corresponding to rotators occur for t esc < 10. We 
will explain this below in the course of the discussion of the HVE. Apart from this narrow, 
serpentine chain, the values of the escape times for F > are scattered mainly over a 
rectangular area with 10 < t esc < 10 4 . This area becomes wider in t with increasing C: 
10 < t esc < 10 3 corresponds to C — 0.01 and 10 < t esc < 10 4 corresponds to C = 0.10. 
These are the particles that are associated with the initial fast decay. The widening of this 
area can be explained with our model from section IIIID1 In the case C = 0.01 (inset of 
Fig. ED the particles with F > escape, similar to the static case with an exponential 
rate, and after a certain time, e.g. t esc ~ 1000, most of them are escaped. Consequently, 
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FIG. 12: Same as Fig. M for the HVE, C = 0.25 and C = 0.05 (inset). 



the algebraic decay establishes itself, see Fig. HI With increasing amplitude, horizontal and 
vertical processes lead to larger changes AF . Due to multiple separatrix-crossing scattering 
the available range in the time t to escape clearly becomes larger. The appearance of rather 
high densities at F ~ 0.9 and F ~ in Fig. [11] and at F ~ 1 in the inset, is explained 
below, in section IIVCI 

For values of F < 0, corresponding to librator orbits, the values of the escape time are 
grouped in an inclined band, i.e. for smaller values of F, the escape time is on average 
higher. In the case C = 0.01, the band stops at F « —0.20, orbits with smaller initial 
values of F just did not escape until 5 • 10 4 collisions were reached. For C = 0.10, this 
band covers almost the hole range in F, due to the larger driving amplitude and the larger 
effective changes AF. 

In Fig. [T21 the results of the HVE are shown. We can match perfectly the exponential 
short-time behavior and the algebraic tail with the two major areas in picture. All particles 
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with initial values F > 0, corresponding to rotator orbits, have escape times of t esc < 10, 
whereas particles on librator orbits with initial values of F < possess escape times t esc > 10. 
Due to high particle velocities, the effects of the horizontal and vertical processes are rather 
small, it takes around t ~ 10 (corresponds to approximately 500 collisions) until the first 
librators are destroyed. Up to 2.5 x 10 4 collisions, only particles with F > escaped 
(C = 0.05), whereas for C = 0.25 particles with F > —0.5 decayed. 

For values F < and t esc > 10 horizontal, narrow layers can be observed in Fig. [121 
The vertical spacing of these layers is 2tt, which is again the period of the breathing ellipse. 
The mechanism at work is the previously (section IIIIBI) mentioned one: When the ellipse 
is expanding, librators are turned into rotators, which can then escape, whereas during the 
contraction period, the rotators are stabilized. Exemplarily, the average escape time as a 
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FIG. 13: Average escape time and standard deviation in the IVE as a function of the initial value 
of F((p ,po) for C = 0.10. 



function of the initial value of F is shown in Fig. [13] for C = 0.10 and the IVE. Starting from 
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F around —2, the escape time is decreasing with increasing F, until the separatrix (F = 0) 
is reached. For values of F bigger than zero, corresponding to rotators, the escape time 
stays approximately constant. Nevertheless, the escape time for a single trajectory with a 
certain initial value F can deviate significantly form the curve shown in Fig. [121 since the 
standard deviation, also shown in Fig. [T31 is quiet large, especially for values of F > 0. 

C. Density Distributions of the Escaping/Nonescaping PAM 



C = 0.01 
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FIG. 14: Density distribution of F((po,po) in the IVE of the whole, the escaped and the remaining 
ensemble for C = 0.01 and C = 0.10 (inset). 

The choice of the initial conditions described in section IIIIBI leads to an amplitude- 
dependent, non-uniform distribution in the density p(F(<pi,pi)) after the first collision 
at (ipi, pi). In Fig. [TU the density distribution of the initial values of F((p,p) for 
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FIG. 15: Same as Fig. [Ljfor the HVE, C = 0.05 and C = 0.25 (inset). 

C = 0.01 and 0.10 is shown (IVE). The three different curves in each of the figures cor- 
respond to the initial (t = 0) [37| density distribution of F(ip,p) of the whole ensemble, to 
the initial density distribution of F for the escaped and to the density distribution of F for 
the remaining (i.e. non-escaped) particles: 

P^max 

Pall(Fo) = Prem(F ) + Pe8C (F ), / P (F )dF = 1. (24) 



J F - 

If we compare these figures with Fig. [TH there is a perfect correspondence between the 
peaks of the initial density distribution of F(ip,p) and the high density regions in the latter 
Fig., i.e. the high density regions are due to the non- uniformity of the distribution of the 
initial values of F. 

With increasing amplitude, the available range of initial values of F still making to an 
escape of the particles possible becomes larger, for C = 0.01 only particles with F > —0.3 
escaped, whereas for C — 0.10 already particles with initial F > —2.4 escaped. 
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The results of the HVE are very similar, see Fig. [12j The main difference is, that 
due to the small effect of a single scattering process, the transition between escaping and 
nonescaping particles in F-space is much sharper and shifted towards higher values of F 
compared to the IVE. Furthermore, the density distribution of the initial values of F explains 
the reverse ordering of No(t = 50) observed in section IIII Bl For C = 0.05, due to our 
definition of the initial ensemble at the innermost ellipse boundary, there are much more 
particles with initial values F ;$ 1, than in the case C = 0.25. Overall, the fraction of 
particles starting on rotator orbits in the case C = 0.05 is larger than in the case 0.25, and 
these rotators will escape fast. 

V. VELOCITY 

In the static ellipse, see section [Til there are two constants of motion. One is the product 
of the angular momenta around the two focus points F(<p,p), which we just studied in the 
context of the driven ellipse in the previous section, the other is the energy. Since the 
potential is constant inside the ellipse, it is sufficient to consider the kinetic energy only. 
The energy of a single particle in the ellipse is given by E to t a i = F^ n = mv 2 /2. Energy 
conservation in the static ellipse thus means \v\ = const. Since all particles have the same 
mass it is sufficient to consider \v\ instead of E to tai- 

When examining F(<p,p), it is instructive to calculate F from the initial conditions 
(ip , po). Doing the same in the case of the velocity \v(t)\ is meaningless, since we know 
|v(0)| = 1 (IVE) or |u(0)| = 100 (HVE) for all particles. Instead, we consider the velocity 
of the particles when they are actually escaping, i.e. |u(i e sc)|, the escape velocity. 

A. Escape Velocity versus Escape Time 

In Fig. [16] distributions of escape times as functions of the escape velocity are shown for 
C = 0.01 and C = 0.10. A point for every pair (\v(t esc )\, t esc ) is plotted in the plane and only 
escaped particles are considered (IVE). At \v\ = 1, there is a vertical line for t < 4. This line 
corresponds to the particles that escape from the billiard without a single boundary collision 
and possess thus an unchanged energy. For times t < 10 2 , the pairs (|u(t esc )|, t esc ) lie on a 
narrow serpentine band. The vertical spacing, i.e. the period of the band is approximately 
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FIG. 16: Escape time versus the absolute value of the escape velocity in the IVE for C = 0.01 and 
C = 010 (inset), each point corresponds to a pair (\v(t esc )\, t esc ). 



2tt, which is the period of the driven ellipse. The band structure is much more pronounced 
in the case of the HVE, see Fig. [T71 where it dominates the overall distribution, naturally 
emanating from \v\ = 100. This correlation between the escape time and the escape velocity 
for small values of t esc can be explained in the following way. The ellipse starts at t — 
from its neutral position with an expanding motion. As long as the ellipse is expanding, 
each time a particle hits the boundary it loses energy and its velocity is reduced. Since the 
particles move very fast compared to the motion of the boundary (HVE), they accumulate 
a lot of collisions until the ellipse reaches its maximal extension and starts contracting. The 
more collisions a particle cumulates during the expansion period, the bigger is the total 
energy loss. The ellipse reaches its turning point at t = 7r/2, i.e. every particle with an 
escape time t esc < 7r/2, will have an escape velocity |u(t esc )| < \v \ = 100. From t = n/2 
on, the corresponding escape velocities will increase until t esc = 3ii/2 is reached, since the 
ellipse is contracting during this time period and every collision with the boundary will 
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FIG. 17: Same as Fig. US for the HVE, C = 0.05. 



increase the energy of the reflected particle. This process is continued until all rotators have 
escaped, which is the case at t ~ 10. This explanation holds for the HVE. Since in the case 
of the IVE, the particle velocities are similar to the boundary velocity, this effect is much 
less pronounced. Nevertheless, it is still visible and mainly due to orbits with initial values 
F m 1, since these orbits skip along the ellipse, accumulating many boundary collisions 
within a short period of time. The main difference of the distributions of the escape time 
for the two above-investigated ensembles is that in the case of the HVE all rotator orbits 
lie on the serpentine band, whereas in the case of the IVE, only rotators far away 
from the separatrix contribute. 

For intermediate times 10 2 < t esc < 10 3 (IVE), the corresponding escape velocities lie 
closely around one for C = 0.01, see Fig. [161 Since the driving amplitude is very small in 
this case, the energy of the particles is not much changed. For larger escape times, the values 
of | v (t esc ) | are a little bit more scattered, since the particles accumulated several boundary 
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collisions resulting in an effective change of |u(i esc )|, still the values deviate no more than 
10% from the the initial value \v\ = 1. As expected, the distribution of the values of |v(t esc )| 
is broadened for larger values of the driving amplitude C, see inset of Fig. [TB] (C = 0.10), 
but the energy gain remains bounded, |u(£ e sc)| ^ 3 even for C = 0.30 (not shown here). 

In the HVE, for escape times higher than 10, corresponding to particles starting originally 
on librator orbits, almost all escape velocities are smaller than the initial velocity \vq\ = 100. 
Horizontal and vertical processes can scatter a particle, moving on a librator orbit, onto 
a rotator orbit (and vice versa), which is a necessary condition for escaping. In which 
direction (towards or away from the elliptic fixed points) a particle is scattered depends 
on many parameters, see section IIII CI but at least the vertical process scatters particles 
during the expansion period of the ellipse always towards rotator orbits. As a consequence, 
particles that turn from librator into rotator orbits accumulate collisions that effectively 
reduce their velocity and scattered them beyond the separatrix, thereby explaining the low 
escape velocities of particles with escape times t esc > 10. One might think of using this 
mechanism to slow down particles. 

B. Distributions of the Escape Velocities 

In Fig. dHl the distribution of the escape velocities is shown for different values of C 
(IVE). With increasing amplitude, the mean escape velocity < |«(t esc )| > is shifted towards 
larger values. Since on average there are slightly more collisions with the contracting ellipse 
than with the expanding one, a larger driving amplitude leads to an increased mean energy 
of the particles. The sharp peak at \v\ = 1 corresponds to the particles that leave the billiard 
without a single collision with the boundary, thus |v esc | = |vo| = 1- 

In the case of the HVE, see Fig. [191 the distribution looks quite differently The majority 
of the particles have an escape velocity smaller than |«o| = 100 and especially at the lowest 
accessible escape velocity \v esc \ « 97 there is a large peak. This is due to the fact that the 
ellipse starts with an expanding motion which deprives the particles energy upon boundary 
collisions and a large fraction of particles decays during that first expansion period. Around 
t = it/2, in the vicinity of the first turning point, the ellipse stays comparatively long 
(boundary velocity ~ 0), more particles escape, leading to the large peak at \v\ ~ 97. 
The asymmetric shape of the distribution is additionally reinforced since the librators that 
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FIG. 18: Distribution of the escape velocity in the IVE for different values of C. 

escape have also low energies. The particles that are not escaped after 2.5 • 10 4 collisions 
accumulated collisions during the expanding and contracting motion of the ellipse and the 
fluctuations in the energy transfer lead to a roughly uniform distribution. The distributions 
for higher values of C look very much like the one shown in the case C = 0.05, except that 
they get wider with increasing amplitude. 

Since most of the particles leaving the ellipse have velocities smaller than \ vq\, the question 
arises, whether the billiard could be used for systematic velocity lowering. To enforce this 
effect, one could try e.g. to choose asymmetric driving laws. We point out that the lowered 
energies of the escaping particles is a feature of the dynamics of the ellipse. In general 
(concerning other geometries) particles are more likely to strike a contracting than a receding 
boundary, which leads one to expect increased energies of the escaping particles. 
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FIG. 19: Distribution of the escape velocity for the HVE, C = 0.05. 



VI. CONCLUDING REMARKS 



We investigated the classical dynamics of the static and especially the driven elliptical 
billiard with an emphasis on the escape rate of an ensemble of particles. As predicted in 
ref. {5 1 in a general context for integrable billiards, we found an algebraic decay in the 
long-time behavior of the static ellipse, due to the integrable structure of the underlying 
dynamics. Besides the energy, the product of the angular momenta (PAM) F((p,p) about 
the two foci is preserved. The sign of the initial value of F determines whether a particle 
moves on a rotator or librator orbit and only the rotators are always (for all hole positions) 
connected with the hole. Consequently, the decay approaches a saturation value N s (e), 
which is maximal for the hole lying at the short side of the ellipse; at this hole position none 
of the librators are connected with it. N s (e) depends on the numerical eccentricity e of the 
ellipse and we predicted this dependence very accurately from theoretical considerations. As 
a consequence, varying e allows us to control the number of emitted particles. 
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When applying harmonic boundary oscillations, neither the energy nor F(ip,p) will re- 
main a constant of the motion. We performed numerical simulations for two different en- 
sembles, representing the two important borderline cases. Firstly, the intermediate velocity 
ensemble (IVE), where |i? | ~ (ouA being the boundary velocity), and secondly the high 
velocity ensemble (HVE), where \vo\ >> ujA. In both cases we observed an initial fast decay 
with an ensuing transition period followed by a non-vanishing (even for large times) near 
algebraic decay. The emission rate depends monotonically on the driving amplitude. The 
changes of F(ip,p) of particles upon a single collision with the boundary are much smaller 
in the case of the HVE, due to the high velocities of the particles. As a consequence the 
resulting decay is similar to the one of the static system. The observed disappearance of the 
saturation value in both ensembles is due to the gradual destruction of the librator orbits 
caused by two fundamental processes: The vertical processes, where upon collisions momen- 
tum normal to the boundary is transferred, making changes in the sign of F(ip,p) possible; 
and the horizontal processes where the particle hits the ellipse due to the boundary motion 
at a different position (compared to the static case) leading again to changes in F({p,p) that 
can result in the transition of a librator into a rotator. Just like in the static system, parti- 
cles starting on rotator orbits (F({p ,p ) < 0) cause the initial fast decay. With increasing 
time, more and more particles with initial conditions closer and closer to the elliptic fixed 
points can escape, due to the just described vertical and horizontal processes, and cause 
the non- vanishing emission rate in the long-time behavior of the decay. We confirmed this 
just displayed strong connection between the escape time and F(ip ,po) by analyzing this 
quantity carefully. In the HVE, the escape rate as well as correlations of the escape time 
and the PAM are modulated with the same period as the ellipse breathes, the ellipse acts 
as a pulsed source. 

Concerning escape velocities, an astonishing feature is observed in the case of the HVE, 
the distribution of the escape velocities is highly asymmetric and particles escape mainly 
with \v(t esc )\ < | v |, the driven ellipse could be used for systematic cooling. To avoid escape 
velocities bigger than |t> |, the use of a point source as an initial ensemble seems reasonable. 
Simulations with thermal ensembles suggested the ellipse as a state transformer, thermal 
ensembles were changed into non-thermal ones. Furthermore, the ellipse could be used as a 
controlable source of particles: if a certain emission rate is is required, this can be achieved 
by tuning the driving amplitude, whereas the numerical eccentricity e of the static ellipse 
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allows us to emit a certain number of particles. 
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